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ABSTRACT 

Context. The F8 star HD 181906 (effective temperature ~ 6300 K) was observed for 156 days by the CoRoT satellite during the first 
long run in the centre direction. Analysis of the data reveals a spectrum of solar-like acoustic oscillations. However, the faintness of 
the target (mv=7.65) means the signal-to-noise (S/N) in the acoustic modes is quite low, and this low S/N leads to complications in 
the analysis. 

Aims. To extract global variables of the star as well as key parameters of the p modes observed in the power spectrum of the lightcurve. 
Methods. The power spectrum of the lightcurve, a wavelet transform and spot fitting have been used to obtain the average rotation rate 
of the star and its inclination angle. Then, the autocorrelation of the power spectrum and the power spectrum of the power spectrum 
were used to properly determine the large separation. Finally, estimations of the mode parameters have been done by maximizing the 
likelihood of a global fit, where several modes were fit simultaneously. 

Results. We have been able to infer the mean surface rotation rate of the star (~4 pHz) with indications of the presence of surface 
differential rotation, the large separation of the p modes (~87 pHz), and therefore also the “ridges” corresponding to overtones of the 
acoustic modes. 

Key words, methods: statistical - methods: observational - stars: oscillations - stars: individual: HD 181906 



1. Introduction 



CoRoT (Convection, Rotation and planetary Transits) is a min- 
isatellite launched on December 26, 2006 and developed by 
the French Space agency CNES (Centre National d’Etudes 
Spatiales) in collaboration with the Science Programmes of ESA 
(European Space Agency), Austria, Belgium, Brazil, Germany 
and Spain. The main objectives are to detect exoplanets and to 



probe the interiors of stars using asteroseismology ( Baglin et al. 
|2006| l. The high-perfomance photometric data offer an unprece- 
dented opportunity to detect extremely low amplitude oscilla- 
tions in many stars ( Michel et al.|2008b) l, in particular, in solar- 
like oscillating stars. The main programme of the CoRoT mis- 
sion consists of 5-month-long runs during which ten designated 
targets are almost continuously monitored with a 32-s cadence 
for seismic studies. Several solar-like oscillation stars have al- 



* CoRoT (Convection, Rotation and planetary Transits) is a min- 
isatellite developed by the French Space agency CNES in collaboration 
with the Science Programs of ESA, Austria, Belgium, Brazil, Germany 
and Spain. 



ready been observed with CoRoT ( |Garcfa et al.||2008] l. One of 
them is HD 181906, a faint F8 star which has been observed dur- 
ing the first long run looking into the galactic centre direction at 
the same time as HD 181420 ( |Barban et al.|200^ . 

In the present paper we report the first detailed seismic anal- 
ysis of HD 181906 using CoRoT data. We first present some 
classical characteristics of the star and infer some expected pul- 
sation properties in Sect. after a description of the CoRoT 
observations (Sect. 3), we derive the surface rotation of the star 
(Sect. 4), we then determine the region where the p modes are 
(Sect. 5) and we explain how to take into account the convective 
background (Sect. 6). We finish by analysing the p-mode spec- 
trum (Sect. 7) and extracting the p-mode parameters (Sect. 8). 

2. Stellar parameters and estimated oscillation 
properties 

The star HD 181906 (or HIP 95221) is known as an F8 dwarf 
with a magnitude m,, = 7.65 and is probably a binary system 
as suggested by astrometric measurements ([Makarov & Kaplan] 
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|2005| [Frankowski et al.|[2007[ l or high-resolution spectroscopy 
([Bruntt 2009|l. The new reduction of Hipparcos data by ' 



Leeuwen 1 2067|l provides a parallax n — 14.72 ± 0.91 mas (i.e. 



a distance c/ = 68 ± 4 pc), which leads to an absolute magnitude 
of My = 3.49 + 0.13. Using BC y ^ -0.044 + 0.058 f or the bolo- 
metric correction in the V band (jMasana et al.|2006ll, we deduce 



the stellar luminosity L - 3.32 + 0.45 Lq. Very recently Bruntt 
( 2009| l has derived global parameters for this star from high- 
resolution spectroscopy. He found a value of Teff = 6300+ 150 K 
which is in perfect agreement with the value he derived from 
2MASS photometry (Teff = 6360 + 100 K). He also deduced a 
surface gravity logg = 4.220 + 0.056 and a metallicity [Fe/H] 
= -0.11+ 0.14. Using his observations as inputs for stellar mod- 
elling, he estimated for this star a mass M - 1.144 + 0.1 19 Mq, 
a radius R - 1.392 + 0.054 Rq and an age of 4.2 +1.6 Gyr. 
These values have been hereafter used to approximate seismic 
quantities and are summarized Table [T] 



Table 1. Global parameters of HD 181906 used for this work. 



n 


14.72 ± 0.91 mas 


L/L© 


3.32 ± 0.45 


Teff 


6300 ± 150 K 


logg 


4.220 ± 0.056 


[Fe/H] 


-0.11 ±0.14 dex 


V sin i 


10+1 kms^' 


M/M© 


1.144 + 0.119 


R/Rq 


1.392 + 0.054 



The frequency of the maximum of the p modes scales as 
the cut-off frequency (jKjeldsen & Bedding 1995 1. This law 



has lately been verified using spectrometric data (Bedding & 



|Kjeldsen|2003|l. Therefore we can express Vmax as follows: 



^ max ~ 



MIMq 



(R/R©)2 ■ 




■ 3050 iMz . 



( 2 ) 



We obtain Vmax ~ 1725 + 225 //Hz 

Finally, the maximum expected amplitude (an estimate of the 
intrinsic mode amplitude in terms of bolometric intensity fluctu- 
ations) can be deduced from|Samadi et al.](|2007b|l: 




( lhq y 




( 3 ) 



where we have changed the solar value of Aq^ux - 2.6 ppm to 
2.53 + 0.11 ppm obtained by calibrating different helioseismic 
measurements of the VIRGO/SoHO package (POM6 & SPM) 
( Michel et al.|2009) l. This gives a value of Amax ~ 5.1 + 0.6 ppm 
for HD 181906. This relation is the combination of two different 



ones: the adiabatic relation proposed by Kjeldsen & Bedding 
( 1995|l to relate mode amplitudes in intensity to mode ampli- 



tudes in velocity and the scaling law proposed by Samadi et al. 
( 2007b | l which gives the mode amplitudes in velocity as a func- 
tion of [LIM). 



We have compared these values to those derived in pre- 
vious works using differents techniques. This star is included 
in the Geneva-Copenhagen survey of the Solar neighbourhood 
( [Nordstrom et al. ( 2004| l revisited by Holmberg et al. (2007 1). 
This survey, mainly based on Stromgren photometry, provides 
an effective temperature Teff = 6382 + 91 K and metallic- 
ity [Fe/H] = -0.18 + 0.10 dex. The age is also estimated to 
2.7/q4 Gyr and the mass to 1.22/qq® Mq. All of these values are 



in very good agreement with those of Bruntt ( 2009)l . By using 
V and 2MASS IR photometry, [Mas ana et al. ( 2006| l recover a 
slightly higher effective temperature of 6532 + 66 K and provide 
an angular semi-diameter SD = 0.102 + 0.001 mas. Thus, us- 
ing the new Hipparcos parallax, we estim ate the stellar ra dius to 



dipp; 

?(£l 



1.50 + 0.10 RfH which is consistent with Bruntt 



(20091. 



For the the rotation velocity, Nordstrom et al. (|2004)l have re- 



ported a measurement based on earlier CORAVEL observations 
and have estimated the rotation velocity vsin/ ^16+1 kms^'. 
However, according to jBruntt ( [2009 ) l, this value is overestimated 
due to the blend with a second star. By taking into account the 
presence of this second star, vsin/ decreases to 10 + 1 kms^'. 
For consistency reasons, we have considered the latter. 

Using these parameters, we are able to obtain approximate 
values of some global seismic parameters using scaling laws. 
The large separation, Av, can be estimated from the mass and 
the radius as follows (Kjeldsen & Bedding|1995|l: 



Wo/ 



R 



■ 3/2 



135 //Hz . 



( 1 ) 



A value of Av 88 + 7 //Hz has been obtained. 

* this number is 5% lower that the one derived by [Masana et al. 
( [2006] , since the parallax has been revised, but is fully compatible 
within the error bar. 



3. CoRoT Observations 



156.6 days of continuous data - from 2007 May 11 until 2007 
October 15 - have been collected with an overall duty cycle 
of 89.3 %. Most of the gaps (each of a few minutes duration) 
are due to data loss during the crossings of the South Atlantic 
Anomaly (see for a detailed explanation: Auvergne et al.|2009 1 . 
These gaps have been linearly interpolated in the light curve to 
avoid putting zeros. We have verified that this interpolation does 
not introduce any spurious frequencies in the Fourier domain. 
The original 1-s cadence raw data have been corrected and cali- 
brated into level-2 (or N2) data following the methods described 
I Samadi et al. (2007a|l. The N2 light curve is sampled on a reg 



ular grid in the heliocentric frame with a cadence of 32s. Then 
we have removed a low-frequency trend due to the aging of the 
CCD ( Auvergne et al.|2009| l, and Anally we have removed some 
outlier points (0.013%). The resultant lightcurve shows a small 
modulation during the first 60 days (see Fig. and then a very 
flat behavior. At this point, it is not possible to disentangle a 
real modulation of the star from an instrumental effect. A sec- 
ond modulation - of around 3 days - is also visible and could be 
related to the surface rotation due to the signature of magnetic 
activity on the surface of the star. 

To compute the power spectrum density (PSD) we have used 
a standard Fast Fourier Transform algorithm and normalized it 
as the so-called one-sided power spectral density (Press et al. 
1992 1. The resultant PSD is shown in Fig.|^ 

In the PSD, several peaks rise above the general trend dom- 
inated by the photon noise from the Nyquist frequency down to 
approximately 100 //Hz, and by the convective noise from this 
frequency until ~10 //Hz ( Michel et al.||2008a) l. Below this, the 
spectrum is dominated by two peaks, at around 4 and 8 //Hz, that 
could be the signature of the surface rotation of the star - already 
seen in its lightcurve - as well as its first harmonic (red dotted 
lines in Fig|^. We will discuss in detail this rotation rate in the 



Garcia et aL: HD 181906 



3 



4000 




-4000 Li . . I I . . 1 . I I I . . I I . i . I I . i 1 . I . 1 . . .L 

0 20 40 60 00 100 120 140 

Time (days) 

Fig. 1. N2-Light curve (in ppm) corrected for the aging of the 
CCD and interpolated onto a regular grid in the heliocentric 
frame. 
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Fig. 2. Power spectrum density of the full 156-day N2-light 
curve shown in Fig[T| The two red-dotted vertical lines show the 
lowest significant peak and the first harmonic that could be the 
signature of the surface rotation of the star. The five blue-dotted 
vertical lines indicate the first harmonics of the orbital period of 
the satellite. 



next section. The signature of the CoRoT orbital periodicity pro- 
duces a peak at 161.7 yuHz together with several harmonics (blue 
dotted lines in Fig.|^. 

The combination of a lower-than-expected signal-to-noise 
ratio of the oscillation amplitudes ( [Michel et al.|200^ with the 
faintness of the target (m^,-7.65) means that the p-mode hump is 
not clearly visible in Fig. but there is a small excess in power 



around ~2000 //Hz. A more sophisticated treatment is necessary 
to clearly unveil the acoustic spectrum of this star. 

4. Surface rotation 

As we have already mentioned in previous sections, the light 
curve of HD 181906 shows a periodic modulation of about 3 
days that produces two peaks in the PSD. It is interesting to anal- 
yse this periodicity in a more detailed way. 

We have started by calculating a time-period diagram using 
wavelets ( Torrence & Compo|1998| l. The advantage of this tech- 
nique relies on the fact that we use a wave, the Morlet wavelet - 
a sine wave modulated by a Gaussian ( Goupillaud et al.|1984| l - 
which has a finite duration and a specific frequency. By chang- 
ing the frequency of this wavelet and sliding it along our time 
series, we calculate the correlation between the wavelet and the 
data. That enables us to compute the Wavelet Power Spectrum 
(see Fig. [^. Most of the power is concentrated along a horizon- 
tal line centered at ~ 2.8 days. This signal appears to be stronger 
during the first half of the run than during the second period. 



Wavelet Power Spectrum 95% Global 




Fig. 3. Left; Wavelet power spectrum for HD 181906 at low 
frequency. Right; Global Wavelet Power Spectrum where the 
dotted-line represents the 95% confidence level. 

With the Global Wavelet Power Spectrum defined as the hor- 
izontal average of the time-period diagram (see Fig.^(right)), 
we observe that most of the power - more than 99% - is con- 
centrated in this main periodicity at 2.81 days. A smaller peak 
is also visible at 1.41 days but containing much less power. This 
method (successfully tested with numerical simulations and with 
solar data from the GOLF instrument ( jMathur et al.||200^ ) al- 
lows us to disentangle between the peak corresponding to the 
main periodicity and that of the harmonic by the simple visual 
inspection of the PSD where 2 peaks of similar characteristics 
stand at 4.04 and 8.2 //Hz. 

A closer inspection of the rotation period, Prot, in the PSD 
reveals that it is composed of a double structure, with a strong 
peak centered at 4 //Hz - 2.9 days (Fig. - and a smaller one 
at 4.45 //Hz (2.6 days). The fact of having these two peaks in- 
stead of just one may suggest the presence of spots at different 
latitudes with a differential surface rotation. The spot modeling 
done by [Mosser et al.j ( 2009a| l states this point explicitly, indi- 
cating two rotation periods associated with a clear gradient of 
the rotational velocity as a function of spot latitude. Another 
explanation, less probable, is the presence of a hot star in the 
background of HD 181906, as suggested by Bruntt ( 2009| l who 
proposed a detailed analysis of the fundamental parameters of 
CoRoT asteroseismic targets based on high-resolution spectro- 
metric measurements. 

If the 156-day time series is divided into two segments of 
78 days, or into three independent 52-day subsets, a slightly dif- 
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Fig. 4. Zoom on the spectram at low frequency. The power spec- 
trum density has been computed using the total length of data. 



ferent low-frequency structure is observed each time, as can be 
seen in Fig. and Fig.[^ 




Table 2. Length of spectra and corresponding rotational periods. 



Spectrum pHz (days) pHz (days) pHz (days) 



Total spectrum 


4.00 


(2.9) 


4.45 


(2.6) 






First 78 days 


4.00 


(2.9) 


4.45 


(2.6) 


4.90 


(2.4) 


Second 78 days 


4.00 


(2.9) 


4.58 


(2.5) 


5.04 


(2.4) 


First 52 days 


4.00 


(2.9) 


4.45 


(2.6) 


4.90 


(2.4) 


Second 52 days 


3.77 


(3.1) 


4.45 


(2.6) 






Third 52 days 


4.00 


(2.9) 






5.12 


(2.3) 



In Fig.[^ the PSD of the first 78 days is plotted in black, and 
the PSD of the last 78 days is plotted in blue. In Fig.|^ the PSD 
of the first 52 days is plotted in black, the next 52 days is plotted 
in red, and the last 52-day segment is plotted in green. Details 
of the prominent peaks that appear in each of these spectra are 
summarized in Table In this table, we listed the frequencies 
(in juHz) and the respective rotational periods (in days) of the 
peaks that appear at low frequency in the Fourier spectrum of the 
full-length series and the smaller subseries that we have consid- 
ered in the analysis. The peaks are located at rotational periods 
in the range from 2.3 to 3.1 days. The different behavior of the 
peaks (in frequency and amplitude) - behavior that depends on 
the observed period - suggests that they could be the signatures 
of differential rotation on the surface of the star, and not the ro- 
tation period of the secondary star in the binary system (which 
would be expected to produce a stable peak in all the considered 
periods). Therefore, HD 181906 shows a slightly larger differ- 
ential rotation than the Sun, but of the same order of magnitude. 
A more detailed analysis of this behavior will be the topic of a 
future work. 

4.1. Constraining the inciination angie 

The angle of inclination of HD 181906 can be inferred from the 
measured value of (vsini)gt,s and the surface rotation that we 
have just derived; 



Fig. 5. Zoom on the spectrum at low frequency. In black the 
Fourier spectrum of the first 78 days of data. In blue the spectrum 
of the last 78 days. 
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Fig. 6. Zoom on the spectrum at low frequency. In black the 
Fourier spectrum of the first 52 days of data. In red the spectrum 
of the 52 intermediate days of data, and in green the spectrum of 
the last 52 days. 



sin i - 



(v sin i)obs 
2nRvrot 



(4) 



Using a rotation frequency of Vrot - 4.0 + 0.15 //Hz and a 
rotation velocity of (v sin i)obs = 10+1 kms“' we obtain an an- 
gle of / = 24 + 3°. We recall this value is obtained by taking into 
account the blend with a second star ( |Bruntt|2009| l. For compar- 
ison, by considering the older value of 16 + 1 km s ', we would 



obtain an angle of 37.5 + 4.5°. From the spot modeling (Mosser 
et al.|200^^ a value of 45 ± 10° has been found. This parame 



ter is extremely important because it affects the amplitude ratios 
of the components of a multiplet of non-radial modes and it is 
strongly connected to the rotational splitting (see, for example, 
[Ballot et al.|2008| and references therein). 



5. Finding the p-mode region 

To find the region where the p modes are centered, we computed 
the Power Spectrum of the Power Spectrum (PSPS) (or the auto- 
correlation of the PSD) in the region determined by the scaling 
laws (see Fig|^and a more detailed explanation in Sect 7.2). 

This revealed a strong peak at half the large separation, Av/2, 
with the location of the peak implying Av ~ 87 //Hz, which 
is close to the upper limit of the value estimated in Section 2, 
using existing, non-seismic data on the star. In order to narrow 
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Fig. 7. Power Spectrum of the Power Spectrum (PSPS) nor- 
malized by its standard deviation in the region between 1400 
and 2100 //Hz. The highest peak corresponds to 6.35 hours, i.e. 
Av=87.5 + 2.6 yuHz. 



down the range more carefully, we then performed the following 
procedure. 

We moved a 300-//Hz-wide window through the frequency 
range of interest, and computed the PSPS at each location. The 
window was shifted in steps of 33 //Hz. We measured in each 
PSPS the height of the peak at Av/2 relative to the local back- 
ground level. The noise distribution in each PSPS will follow 
two-degrees-of-freedom statistics; given this known distribution, 
it is possible to calculate a false-alarm probability for the Av/2 
peak to appear by chance (as part of the background) anywhere 
in the PSPS. We calculated the level corresponding to a 5% of 
probability of appearing by chance ( [Chaplin et al.|200^ . 

The maximum relative height of the Av/2 peak (i.e., the 
height divided by the measured background in the PSPS) is plot- 
ted in Fig. ^ We find here that the p-mode power is prominent 
between 1400 and 2100 //Hz. 




Fig. 8. Maximum of the PSPS in a region around the ex- 
pected half-large separation computed in 300 //Hz-wide win- 
dows shifted every 33 //Hz. The horizontal error bars show the 
range in frequency in the PS that was analyzed (300 //Hz). The 
dashed line represents the 95% conhdence level. 



6. Background fitting 

In order to estimate the non-p-mode background, we fit the fol- 
lowing three-component model to those parts of the PSD where 
the observed p-mode power is insignihcant (PSDbg): 



2 



PSDbg(v) = 



St 



A; 



V !=1 



+ (vB,)C. 



+ D 



(5) 



There are two power-law components in the summation: a 
component to represent the significant power at very low fre- 
quencies from surface activity; and a component to represent 
power from surface granulation. Both power-law components 
are modelled in terms of three parameters; a power spectral den- 
sity, A; a characteristic timescale, B\ and a power-law index, C. 
The third component in Equation D, models the contribution 
from shot noise. 

Rather than fit the raw power spectrum, we fitted a smoothed 
spectrum generated by applying an A^s-bin-wide boxcar, taking 
the independent averages only. With ^ 30, independent N^- 
bin averages show normally-distributed scatter about the (under- 
lying) limit spectrum we seek to estimate from the fit. A standard 
least-squares fitting was therefore applied, with weights fixed by 
the uncertainties on each independent As-bin average. These un- 
certainties were each given by sj y/Nl, where s is the standard 
deviation of the Ns contributing power values. 

The alternative approach is to ht the raw spectrum by max- 
imizing a likelihood function commensurate with the two- 
degrees-of-freedom statistics. It turns out that the smoothed 
spectrum may also be fitted by applying the same maximum like- 
lihood estimator (as was shown by ( |Appourchaux|20()4) l). These 
approaches are, however, more sensitive to the choice of initial 
first-guess parameters than is the least-squares htting approach 
applied here, and can as such be prone to poor convergence. 

There are a total of seven free parameters defining the back- 
ground model in Equation]^ We did not, however, fit them all 
simultanously: some parameters were fixed, as we now go on to 
explain. 

We assumed power in the very low-frequency activity com- 
ponent arises predominantly from the exponential decay of ac- 
tive regions and plage. This decay implies a limiting power den- 
sity spectrum that is Lorentzian, hence we fixed the activity- 
component index to a value Ci = 2 during fitting. 

Exercises performed with the asteroELAG artificial astero- 
seismology data ( [Chaplin et al.|2008 1 showed that attempts to fit 
Al and Bi simultaneously with the other parameters could lead 
to instability in the procedure and poor convergence. We found 
that we could stabilize the fitting by fixing Ai at the value of the 
power spectral density of the first element of the averaged spec- 
trum, leaving Z?i as the only parameter of the activity component 
to be htted. This approach does mean that some care is needed 
in interpreting the best-fitting value of Bi, since the value for 
Al can be affected by the appearence in the power spectrum of 
narrow-band features arising from rotational modulation, which 
are not modelled in Equation]^ The aforementioned approach to 
the htting does however give a good representation of the power 
from the active-region decay that leaks into the frequency region 
where the granulation is important; which in turn means we can 
in principle have more conhdence in the accuracy of the best- 
htting granulation parameters. 



7. The p-mode spectrum 

Eigure plots the PSD of the full-length lightcurve smoothed 
with a boxcar of 70 points (~ 5.2 //Hz). A clear excess of power 



6 



Garci'a et al.: HD 181906 





I 



|:I.l 
^ 1-0 
:0.9 
:O.0 
- 0.7 

ro.6 



0.5 




2000 

Frequency (a^Hz) 



Fig. 9. Power spectrum density of the full length light curve 
smoothed by a boxcar of 70 points (Top). The continuous green 
curve is the background fitted using the method explained in 
Section]^ The bottom panel shows a closer region of the p-mode 
band between 1000 and 3000 fjHz after subtracting the back- 
ground (green curve in the top panel). The blue curve in both 
panels is the PSD but smoothed using a boxcar of 5000 points. 



is observed, relative to the best-fitting background model (green 
curve), in the region where the p-mode excess was detected (see 
Sect. 1^. To visually enhance the excess we also show in blue a 
heavily smoothed spectrum computed by applying a 5000-point 
boxcar (~ 370 juHz) to the raw spectrum. 



Figure [T0| shows the Echelle diagram ( Grec et aL]|1983) l of 
the 70-point-boxcar smoothed PSD. The diagram covers the fre- 
quency region from 792.5 to 3000 fiHz, with each horizontal 
strip covering 87.5 yuHz of the spectrum (see Fig.[7|i. Two ridges 
appear in the diagram, corresponding to the odd and even modes. 
Inspection of the diagram shows that it is not possible to discrim- 
inate visually between the two ridges, hence the angular-degree 
tagging is uncertain. 



7.1. Global amplitude 

The procedure for measuring the mode amplitudes begins by 
averaging the power spectral density in independent frequency 
slices of qAv, where q — 1 or 2, and Av is the estimated large fre- 
quency spacing of the acoustic mode spectrum. The fitted back- 
ground is then subtracted. The resulting residuals in power spec- 



Fig. 10 . Echelle diagram of the full PSD smoothed by a 70 points 
boxcar starting at a frequency of 792.5 pFIz and using a folding 
frequency of 87.5 pWz. 



tral density will be greater than zero over ranges occupied by 
significant mode power. 

Re-calibration of the residual averages then allows for an es- 



timate of the equivalent / = 0 mode amplitudes (Kjeldsen et al. 
|2008[ ). To re-calibrate, one must: (i) multiply by qAv, to give the 
average power across frequency intervals of this length; and (ii) 
normalize to power per I - 0 mode, by dividing by qViot, where 
Vtot is the combined visibility of the I - 0, 1, 2 and 3 modes 
(again, see |Kjeldsen et al.| ( 2008| l). 

In practice we ran a boxcar of width qAv through the spec- 
trum in order to better estimate the maximum power (one could 
miss the true maximum if the independent slices cut through 
pairs of modes). The measured maximum mode rms amplitude 
was Amax = 2.9 + 0.3 ppm. 

We also applied a different approach to estimation of the am- 
plitudes. This approach assumed the envelope of excess power 
due to the modes could be modelled by a Gaussian function. We 
fitted a Gaussian profile to the re-calibrated residual averages 
(see above), and from the best-fitting maximum we estimated 
^max = 2.8 + 0.1 ppm. The location in frequency of maximum 
also allowed us to estimate the frequency of maximum power, 
which was v^ax = 1912 + 47 yuHz. 

To convert these instrumental values into maximum intrinsic 
bolometric amplitudes per radial mode, Aboi(f = 0), we took into 
account the instrumental response functions for CoRoT, as pre- 
sented by ( [Michel et al.||2009 1 . We then found that Aboi(f = 0) 
= 3.26 ± 0.42 ppm (Michel et al. 2008a i. This amplitude is 
about one third smaller than the expected value deduced from 
the scaling laws. Such a deficit is also seen for the other hot, 
solar-like oscillators observed by CoRoT ( Michel et al.|200^ . 
HD 181906 is not significantly undermetallic, so it is hard to 
appeal to a metallicity effect as a possible explanation (e.g., as 
in jMosser et al.|p0081 l). However, we note that this star - like 
the other solar-like CoRoT targets HD 49933, HD 175726, and 
HD 181420 - rotates significantly faster than the Sun (10 times 
for this target) and exhibits clear magnetic activity (cf. Sect. 4). 
It could be an indication that the magnetic activity plays a role 
in suppressing the p-mode amplitudes by changing the surface 
convection and thus modifying the efficiency of p-mode exci- 
tation processes. Some 3D compressible radiative magnetohy- 
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drodynamics simulations Pacoutot et al.|2008] l have shown such 
effects. 

There is another possibility, which is that the discrepancy 
could be explained by a blending effect. Bruntt| ( |2009 1 have 
shown evidence in the spectrum of HD 181906 for the super- 
imposition of a second spectrum, which could be a companion 
or a star in the background field. If this is the case, the relative 
amplitudes of modes could easily be underestimated by a factor 
of 1.5, making them compatible with the observations. 

7.2. The mean large separation 

Different methods have been used for extracting the large sep- 
aration of the p-mode spectrum. A good signature of this large 
separation can be derived from the PSPS, or from the autocorre- 
lation of the power spectrum, computed in the region where the 
excess of power has been detected (see Fig.|^. 

Taking the frequency interval where the p-mode excess has 
been found (see Sect. |^, i.e. from 1400 to 2100 pHz, we find 
that the PSPS is dominated by a peak at 6.35 hours (see Fig. [7]l, 
which corresponds to half the large separation of 43.75 pHz. 
Indeed, the main periodicity we found is not the large separa- 
tion itself but the distance between peaks of even and odd de- 
gree. Therefore, the large separation is Av = 87.5 + 2.6 //Hz (see 
Fig. 0- Depending on the range in frequency used to look for 
the large separation, different groups have found slightly differ- 
ent values. For example, if we reduce the search range to 1400 
to 2000 //Hz, the large separation is reduced to a value of 85.7 + 
2.6 //Hz. 

Another powerful way to estimate the value of the large sep- 
aration is to build the so-called collapsogram ( |Korzennik|2008| l . 
The first step consists in calculating all the Echelle diagrams for 
a frequency range in which we expect the large separation (e.g. 
60-100 //Hz). The second step then consists in collapsing the 
vertical dimension of each Echelle diagram. In the last step the 
collapsed diagrams are stacked vertically, to yield the final col- 
lapsogram. 
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Fig. 11. Collapsogram of the smoothed power spectrum (11 
bins) obtained by collapsing the power spectrum in the fre- 
quency range between 1280 and 2300 //Hz. The vertical axis 
corresponds to the large separation that was explored, i.e., from 
75 to 100 //Hz. The horizontal axis indicates the frequency, from 
0 to 100 //Hz, which is the maximum large separation that was 
looked for. Note that there are two prominent excesses of power 
in the collapsogram, corresponding to a large separation of about 
85 to 90 //Hz, at ^21 and ^^67 //Hz. 



analyzed in this pap er) and HD 181420 (|Barban et al.|2009 1, 60 



days for HD 49933 ( Appourchaux et al.|2008| l and only 27 days 



for HD 175726 ( |Mosser et al.||2009b| l. The data from Procyon 
correspond to velocity measurements obtained during 26 days 
for a multisite campaign involving 1 1 telescopes (|Arentoft et al. 
[ 20081 ). 

The data in Table are obtained mainly from the above 
mentioned papers, except the bolometric maximum amplitude 
per radial mode (A^ax) for the 3 first mentioned CoRoT stars 



us to explore on a single diagram a wide range of possible val- that were obtained from Michel et al. 


P008a[). Eor Procyon, the 


ues of the large separation. Secondlv, it increases the visibilitv of Teff, v sin i, and [Ee/H] were obtainec 


from 


Allende Prieto et al. 



the separation compared to that seen in the Echelle diagram (al- 
though it does rely on the separation being almost uniform over 
the frequency range used to construct the collapsogram). 

In Eig. 



is a translation from cm/s to ppm done in Arentoft et al. ( 2008)l 



by the authors, a value that agrees well with the one obtained by 
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we present the collapsogram of the smoothed 
power spectrum obtained with a weighted moving average over 
1 1 bins. The presence of the two mode ridges is seen, with the 
large separation indicated as being between 85 to 90 //Hz. Note 
that we have considered only the power spectrum in the fre- 
quency range between 1300 and 2300 //Hz to build the collapso- 
gram. 

7.3. Comparison with other observed F stars 

A comparison of published global parameters (and their asso- 
ciated error bars when available) from different E stars is pre- 
sented in Table j^where Av is the large spacing, v„,ax is the fre- 
quency in the spectrum with the maximum power and A„ax is 
the bolometric maximum amplitude per radial mode. In the case 
of Procyon, the bolometric amplitude is obtained from the ve- 
locity measurements following [Arentoft et al.j 



Bruntt et al. (2005 1 using white-light photometric data from the 



[). The stars 

HD 181906, HD 49933, HD 181420 and HD 175726 are ob- 
served by CoRoT. However, the length of observations is not the 
same for all of them, being 156 days for HD 181906 (the star 



WIRE satellite. 



8. Extracting individual p-mode characteristics 

Although fitting low-degree p-mode profiles in helioseismology 
and asteroseismology might appear very similar, the unknown 
stellar inclination angle makes the fitting of asteroseismic data 
much more difficult (see for example Appourchaux et alT[ (2008 1 
for the case of the star HD 49933 or [Gizon & Solania[ ( 2003 i 
for Monte-Carlo simulations with artificial p-mode profiles). It 
is not only the lower signal-to-noise ratio of the p-mode astero- 
seismic signal which makes the fitting difhcult but also the high 
correlation between the inclination and the rotational splitting 
( [Ballot et al.|2006[|2008] l. Because of that, the determination of 
these two parameters can be rather poor and will consequently 
affect the determination of the other parameters (frequencies, 
widths, heights...). Therefore, instead of fitting individually each 
multiplet or pair of modes - as it is commonly done in helio- 
seismology - we choose to perform a global fitting of all the 
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Table 3. Comparison of published global parameters for different F stars. 



Stars 


HD181906 
This paper 


HD49933 

Appourchaux et al. 2008 


HD181420 
Barban et al. 2009 


HD 175726 
Mosser et al. 2009b 


Procyon 

Arentoft et al. 2008 


Spectral Type 


F8 


F5 


F2 


F9/G0 


F5 


Teff 


6300 ± 150 K 


6780 ± 130 K 


6580 ± 105 K 


6000 ± 100 K 


6514 ± 27 K 


[Fe/H] 


-0.11 ±0.14 dex 


-0.37 dex 


0.00 ± 0.06 dex 


-0.22± 0.1 dex 


-0.05 dex 


V sin i 


10 ± 1 kms“* 


9.5- 10.9 km s-‘ 


18 ± 1 kms-' 


13.5 ±0.5 km s-‘ 


3.16 ± 0.5 kms-‘ 


Av 


87.5 ± 2.6 //Hz 


85.9 ±0.15 //Hz 


~ 75 //Hz 


~ 97/zHz 


~ 55 //Hz 


^max 


1900 //Hz 


1760 //Hz 


1500 //Hz 


2000 //Hz 


900 //Hz 


A 

^max 


3.26 ± 0.42 ppm 


4.02 ± 0.57 ppm 


3.82 ± 0.40 ppm 


~ 1.7 ppm 


~ 8.5 ppm 



multiplets above a given amplitude threshold around the maxi- 
mum of the p-mode hump, assuming that the rotational splitting 
is independent of the frequency (see Appourchaux et aL]|2008| 
for all the details). This type of global method was pioneered 
by |Roca Cortes et~aL ( 1999[ ) using solar data. By doing so, the 
splitting and the inclination angle are better constrained, even 
though HD 181906 is then modeled as a rigidly rotating star. 
Each multiplet is described by five parameters: the central fre- 
quencies of the modes I - 0, 1,2, one line width (the same for 
all modes within a large separation), and one mode height. We 
assumed the same visibility ratio between angular degrees as the 
ones used in full-disk helioseismology and the visibilities be- 
tween m-components given by |Gizon & Solanki| ( |2003| ). 

Eight teams performed a global fitting of the HD 181906 
data, seven of them using a maximum likelihood minimization 
and one using a least square fitting over an averaged spectrum. 
All the teams did the fits for both scenarios (A and B) depending 
on the identification of the ridges as the odd or even modes. 

Two teams decided to split the observations in four subseries 
and computed the Joint Power Statistics (JPS, see Appendix A 
for a detailed explanation). The JPS is an alternative method that 
contains the same underlying information as the average spec- 
trum but with a different treatment of the noise. The rest of the 
teams worked on the full resolution power spectrum. The num- 
ber of overtones fitted was 5, 7, 9 and 16 overtones. Different 
strategies have also been followed by each team to obtain the 
initial parameters for the fit. In particular, one team used the re- 
sults from HD 49933 ( [Appourchaux et ak] 2008 1 as a guideline 
since their PSDs look very similar and HD 49933 has a much 
better signal-to-noise ratio. Eigure shows the superposition 
of the PSD of HD 181906 and the one of HD 49933 - prop- 
erly scaled in amplitude - and shifted by 17 //Hz. Thus, the 
initial guesses were obtained directly from the fitted values of 
HD 49933. Another team repeated the fits 200 times for each 
scenario, adding scatter to the initial parameters in order to test 
the robustness of the fits and identify any correlation between 
the initial seed values and the fitted ones. As expected, the in- 
clination angle was identified to be strongly correlated with the 
initial seed angle. 

The mode parameters were then extracted using different ap- 
proaches based on the recipe developed by the Data Analysis 
Team (DAT) and explained in [Appourchaux et ~aL] ( [2008[ l. 
However, in order to have a more reliable characterization of 
the p-mode signal, given the low signal-to-noise ratio of the ob- 
servations, some more a-priori information needed to be intro- 
duced. In order to stabilize the fits and avoid systematic outliers, 
some of the fitting teams fixed some parameters to be equal over 
the range of considered overtones. Eor instance, without some 
a-priori conditions, the minimization procedures failed to re- 
turn reliable estimates at high frequency, fitting mostly spikes 
and over-estimated heights. Tables]^ and show the fitted mode 
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Fig. 12. Smoothed power spectrum of HD 181906 (black curve) 
and smoothed and shifted by 17 //Hz power spectrum of 
HD 49933 (blue curve). 



frequencies obtained from the less constrained fit for both sce- 
narios A and B respectively, the analysis being performed on 
the power spectrum of the whole HD 181906 time series. The 
only condition applied to those fits is that the mode widths and 
mode heights were set to be uniform (the same for all the modes) 
over the fitted frequency range. In the central part of the p-mode 
hump, those frequency estimates are consistent within 3cr with 
the estimates returned using more constrained fits and different 
dataset lengths. In the low- and high-frequency ranges, the dis- 
crepancies are larger but still within 5cr. Unfortunately, neither 
of the scenarios seems to be favoured: for instance, the likeli- 
hood ratio test does not allow us to disentangle between the two 
possibilities, the likelihood ratio being not greater than 5 (for the 
fits presented here). 



Table 4. HD 181906 frequencies (in //Hz) for scenario A. 







Scenario A 




# 


Z = 0 


1 


1 = 2 


1 


1524.80 ± 0.86 


1565.61 ± 0.38 


1521.98 ±0.42 


2 


1604.89 ± 0.25 


1645.98 ± 0.41 


1596.46 ±0.51 


3 


1684.62 ± 0.40 


1734.98 ± 0.53 


1687.32 ±0.51 


4 


1773.14 ±0.27 


1814.29 ±0.31 


1772.43 ± 0.39 


5 


1862.25 ± 0.53 


1898.77 ± 0.27 


1861.74 ±0.39 


6 


1947.42 ± 0.30 


1988.32 ±0.37 


1946.26 ± 0.28 


7 


2035.88 ± 0.62 


2076.45 ±0.51 


2036.30 ± 0.95 
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Table 5. HD 181906 frequencies (in /iHz) for scenario B. 



Scenario B 


# 


1 = 0 


/= 1 


1 = 2 


1 


1565.96 + 0.38 


1521.94 + 0.41 


1560.55 ± 0.65 


2 


1646.63 ± 0.35 


1604.76 ± 0.29 


1644.97 ± 0.41 


3 


1734.59 ± 0.26 


1684.74 ± 0.62 


1735.54 ± 0.53 


4 


1814.83 ± 0.41 


1772.68 ± 0.27 


1814.16 + 0.41 


5 


1898.79 ± 0.43 


1861.91+0.43 


1898.46 ± 0.40 


6 


1988.50 + 0.33 


1947.04 ± 0.25 


1988.72 ± 0.44 


7 


2075.91 ± 0.32 


2036.48 ± 0.69 


2077.20 ±0.38 



Nevertheless, we can still obtain some global characteristics 
of HD 181906. On one hand, the htted inclination angle, for both 
scenarios (around 47° + 3.5°) , is very close to the one derived 
from the analysis of the low-frequency peak of the PSD (see 
Sect. 4) when a (v sin i)obs ~ 16+ 1 km s“' is considered but quite 
far from the ~ 24° obtained considering (vsini)o*i ~ 10 kms^'. 
On the other hand, the extracted rotational splitting (5.8 and 6.1 
+ 0.14 fiHz respectively for A and B scenarios) is overestimated 
compared to the surface rotation rate (see Sect. 4). This could be 
due to an increase of the rotation rate in the stellar interior or a 
biased estimation due to the a-priori conditions applied to the fit 
(a common height and width for all the htted modes).These re- 
sults show, once again, the strong correlation between these two 
parameters and the necessity to have an external good determi- 
nation of the inclination angle or the rotation of the star to be 
able to disentangle between the possible solutions. 

Another parameter - the global line width - was found to 
be rather small (1.25 -0.3/+0.4 and 1.28 -04/+0.6 pHz for sce- 
narios A and B respectively), which is probably a consequence 
of the low SNR that could bias this estimate when we apply 
our peak-bagging codes: the low SNR means the modes will 
appear “spikier” than they really are (we lose the Lorentzian 
tails in the background noise), meaning we will be biased to ht- 
ting small line widths. We expect that the autocorrelation of the 
power spectrum will have better SNR than that seen in individual 
modes and hence we look there for the features associated with 
mode line width. The features in it are quite wide, implying that 
the line widths are large and probably in line with what we saw 
for HD 49933 (see Fig.[T2|). Indeed, the widths of the modes in 
both stars, around the maximum of the p-mode spectrum, seem 
to be similar which favours shorter mode lifetimes than in the 
Sun. 

Finally, we would like to emphasize that the combination of 
a small SNR, a possible small inclination angle of the star (close 
to 25°) and large mode line widths produce similar Lorentzian 
prohles for the modes 1=0, 2 and 1 ( [Ballot et al. 2004 1. This 
effect could contribute to explaining why we obtain very similar 
htting results in both scenarios. More work will be necessary to 
obtain more reliable individual p-mode parameters. 



9. Conclusion and perspectives 

In this paper we have shown the hrst seismic analysis of 
HD 181906 (HIP 95221), a faint CoRoT target (m,,=7.65) that 
has been observed continuously during 156.6 days in 2007. The 
surface rotation of the star has been inferred by analyzing the 
very low-frequency part of the power spectrum. A rotation rate 
of 2.9 days (4 pHz) has been established. The presence of a sec- 
ond peak close to the previous one in the power spectrum and 
produced during a different period of time, has been interpreted 
as the signature of the presence of magnetic structures on the 



stellar surface at different latitudes and, therefore, it has been 
deduced that this star could have a slightly higher differential ro- 
tation than the Sun. Coupling this rotation rate with the previous 
result of the (v sin i)ohs ~ 16 + 1 km s“' has allowed us to infer an 
inclination angle of 37.5° + 4.5° but this value could be reduced 
to 24° + 3° if we consider that HD 181906 has a companion 
which implies a reduction in the (vsinOo^^ to 10+1 kms^'. 

A comb-like structure has been unveiled between 1400 and 
2100 yuHz which corresponds to the acoustic -mode spectrum 
with a maximum power at 1912 + 47 pHz and a large separa- 
tion of 87.5 + 2.6 fjHz measured inside this frequency interval. 
However, the low signal-to-noise ratio of the modes prevents us 
from unambiguously identifying them. To go further, more a- 
priori information is needed to constrain the fits. 
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Appendix A: Joint Power Statistic 



In cases where the signal-to-noise ratio of a power spectrum 
is low, it can be useful to increase the signal and decrease the 
noise at a cost of decreasing the resolution in frequency. This 
has often been accomplished by calculating either the arithmetic 
or geometric mean of a set of power spectra obtained from in- 
dependent contiguous subsets of the complete time series. The 
expected statistics of these averaged power spectra can be calcu- 
lated, however they do not have the same negative exponential 
distribution as a single power spectrum. The standard maximum 
likelihood estimation (MLE) fitting techniques are designed for 
power spectra with negative exponential statistics (i.e. with 
2 d.o.f.), it is therefore desirable to calculate an averaged power 
spectrum with the same statistics as the original spectrum. 

The joint power statistic (IPS) ( Sturrock et al.|2005l l esem- 
bles a correlation function in the sense that it increases the con- 
tribution of signals present in a number of independent power 
spectra while decreasing the uncorrelated noise. It also has the 
important property that the IPS is distributed with negative expo- 
nential statistics and can, therefore, be immediately htted using 
established MLE techniques. The IPS can be calculated using 
the power spectra of any number of independent subseries, with 
the resolution in frequency decreasing by a proportional amount 
compared with the power spectrum of the complete time series. 
In the case of HD 181906, the best compromise between increas- 
ing the amplitude-to-background ratio and maintaining sufficient 
resolution in frequency was found when the complete detrended 
time series was divided into four contiguous subseries. The ap- 
propriate fourth order IPS can be approximated by 



- 



3.881X2 

1.269 -hX’ 



(A.l) 



where X is the geometric mean of the power spectra (5,) calcu- 
lated from four independent subseries 






(A.2) 



In the IPS calculated from four subseries, after smooth- 
ing with a Gaussian filter of width lcr=2yuHz, the average 
amplitude-to-background ratio for mode peaks between 1700 
and 2000 juHz is 2.8. This compares with the equivalent full se- 
ries smoothed PET amplitude-to-background ratio of 1.7 over 
the same range. While the IPS can be htted assuming nega- 
tive exponential statistics, the mode amplitudes obtained will be 
higher than those obtained from a single power spectrum due to 
the increase in power of correlated signals in the IPS. Similarly, 
the stellar background parameters obtained from a IPS ht will be 
lower than those obtained from a power spectrum ht. The mode 
frequencies and linewidths obtained from the IPS are the only 



